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Abstract 

We give out the time evolution solution of simultaneous amplitude and phase damping for any continuous 
variable state. For the simultaneous amplitude and phase damping of a wide class of two- mode entangled 
Gaussian states, two analytical conditions of the separability are given. One is the sufficient condition of 
separability. The other is the condition of PPT separability where the Peres-Horodecki criterion is applied. 
Between the two conditions there may exist bound entanglement. The simplest example is the simultaneous 
amplitude and phase damping of a two-mode squeezed vacuum state. The damped state is non-Gaussian. 

PACS 03.65.Ud 03.65.Yz 03.67.Mn 42.50.Dv 

1 Introduction 

Quantum entanglement or inseparability plays a major role in all branches of quantum information and quantum 
computation. Peres P proposed a criterion for checking the inseparability of a state by introducing the partial 
transpose operation. This condition is necessary and sufficient for some lower dimensional discrete bipartite 
systems but is no longer sufficient for higher dimensions 0. Despite many studies on the discrete states, much 
attentions have been paid to the continuous variable states Recently, quantum teleportation of coherent 
states has been experimentally realized by exploiting a two-mode squeezed vacuum state as an entanglement 
resource^]. Due to the decohence of the environment, a pure entanglement state will become mixed. Thus it 
is important to know if a given bipartite continuous variable state is entangled or not. The decoherence may 
be caused by coupling to the thermal noise of the environment, amplitude damping, quantum dissipation and 
phase damping. Besides the phase damping, the other three types of decoherence preserve gaussian property 
of the state, and a two-mode squeezed vacuum state will evolve to a two mode gaussian mixed state. For the 
separability of two mode gaussian state, the positivity of the partially transposed state is necessary and sufficient 
[H]|E]|Z|' However, a gaussian state will evolve to a non-gaussian one by phase damping, and the case of a two 
mode squeezed vacuum state under the only decoherence of phase damping was perfectly solved(Sj- In real 
experiments, the general situation which should be taken into account is the coexistence of noise, amplitude and 
phase damping. Theoretically, the separability and entanglement of non-gaussian state are seldom investigated, 
here we provide an example. 

2 Time Evolution of Characteristic Function 

Considering the simultaneous damping, the density matrix obeys the following master equation in the interaction 
picture ^ = [Ci+C2)p. Where Ci is the amplitude damping part 

Cip = ^ yK"'^ + l)(2aipa+ - a^Uip- pa^ai) +ni[2af pa.i - aiO^ p - paiol)], (1) 

i 

with n the average photon number of the thermal environment. And C2 is the phase damping part (e.g. 
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The state can be equivalently specified by its characteristic function. Every operator A £ B{Tl) is completely 
determined by its characteristic function XA ■= tr[A'D{^)] [H], where = exp(^a+ — /x*a) is the displacement 
operator, with pi — [pti, pi2, ■ ■ ■ , Hs]'^ : a = [ai, ^2, • • • , ^s] and the total number of modes is s. It follows that A 
may be written in terms of XA as ll^: A = /[Hi ^Sr^]X-4(M)^(~A*)- The density matrix p can be expressed 
with its characteristic function x- The amplitude damping equation of x and its solution are well known, they 
are ^ = -^J2i'^ii\f^t\ + Im^I^x), xifJ-^t) = xifJ-i^'^ ,0) exp[-J2^in^ + ^)(1 - e-^.t) We now give 

out the phase damping equation of it will be 

if we denote fii as \fii \ e^^\ We can see that with the characteristic function the amplitude damping equation is 
described by the amplitude of the parameter /i^, the phase damping equation is described by the phase of the 
parameter pLi. The solution to the phase equation of x then will be x (Mi I^* t^) = Yii (27r7ii)~^^^ J exp(— 
X (yLte", ^*e~", O) dx, where /ze" is the abbreviation of [/Ltie*^^ , ^2e"^, • ■ ■ , /Xse"=].The simultaneous amplitude 
and phase damping to any initial characteristic function then will be 

X (Ai, t) = n i^^^ty'^' J exp Et-^ + ("^ + - e""^'*) l/^^l']^ , /.^e-^-^N o) dx. (4) 

i i 

The density matrix then can be obtained as well by making use of operator integral. 

We will concentrate on the simultaneous amplitude and phase damping of two-mode x—p symmetric Gaussain 
stateJU whose characteristic function is x(/^,0) = exp[— (Aio |mi|^ + ^20 \^^2^) + BqPi^12 + B^p\p2]- ■ The time 
evolution of x is 

x(m,M*,0- / exp[-K|/ii|VA^|/i2|') + i?MiM2e" + i3*MlM;e-")]. (5) 

Where A[ = e^^^^'A^o + (n^ + ~ e"'^'*), B = e^^^^^+^'^'^^B^,^ = ^(71 +72). The integral on x is one fold, for 

simplicity, we denote ^^^—^ J S:xp[—-^]f{x)dx as f{x). The x—p symmetric Gaussain state set is a quite large 

set. It contains two-mode squeezed vacuum state |^) = j^j^^hT X]n(^^^^'')" I""' ( is the squeezing parameter) 
and two- mode squeezed thermal state [2] as its special cases, with Aio = A20 ~ \ cosh2r, Bq — ^ sinh2r and 
^10 = A20 = {no + ^) cosh2r, Bq ~ (no -I- |) sinh2r respectively. 



3 PPT separability 

For the sake of simplicity in description, let us firstly consider the situation of A[ = A'2 = A' , B = B* . The 
general case will be obtained straightforward and be described at the end of this section. The first question is 
that if the state after damping is entangled or not. Then how much is the entanglement left? The necessary 
condition of a bipartite state being entangled is that the partial transpose of the density operator is not positive 
definite 0. The partial transpose operation changes the characteristic function in the fashion of : x M2) =^ 
x(/ui,— /X2) = x^"^ (a*1) M2) • For the separability of a non-gaussian bipartite state, a necessary condition was 
proposed by Simon|S] in terms of the second moment of the state. In the original literature canonical operators 
were used, here we use creation and annihilation operators instead. The necessary condition comes from the 
non-negativity of p^"^ and the commutation relations. For any Q — rjrf^ with 77 = ciOi -I- C2CL2 + c^cli + 040^ 
of every set of complex coefficients q, one has (Q) = tr{Q p^^) > 0, Hence the second moment matrix of p^^ 
should be semi-positive definite. The second moment such as tr{p^'^ af aj) can be obtained from the second 
derivative of x with respect to pi and p* , we have 

A-l> Be-^^ (6) 

Where A = A' + ^ . Other necessary conditions may come from when rj is the linear combination of higher 
power of the creation and annihilation operators, and they may be tighter than Ineq.®. And this is really the 
case. 
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We will find a tighter condition by exploring the negative eigenvalues of the partial transpose of the den- 
sity operator. The eigenequation of p^^ can be simplified as the eigenequations of a serials of matrices (see 
Appendix) . 

- - ^^)-^^"E (r;) (:) (^) 

where C = 1 — -^tztb^j D — j^-i^^i ■ When m = 0, one has the first eigenvalue A^°' — {A^ — B^)^^ which is 
always positive. The matrix M^™) possesses the symmetry of Af;*-™-* = M^^^ , so that it can be reduced, 
and we get more analytical solutions. The negative eigenvalues may appear at Aq^^* — (A^ — B'^)^^{C — De~'^*), 
Ag^'' = (A"^ — B^)^^{C^ — D'^e~'^'^*) and so on. Hence one of the necessary conditions of the non- negativity of 
p^'^{t), so that the necessary condition of a damped state p{t) is PPT separable is that 

C>De-''K (8) 

We will prove that this condition is also sufficient for PPT separability. We turn to the detail properties of 
matrix MW. The necessary condition of separability comes from M(2) is C > De'^'^, this is a trivial result 
compared with Ineq. (jHJ We have checked other solvable eigenvalues for necessary condition of separability. They 
are also trivial compared with Ineq.®. It may be anticipated that separable conditions come from all other M'™^ 
are weaker than Ineq.jHl, that is Ineq.ljHJ is also a sufficient condition of PPT separability. To prove this, we just 
need to consider the PPT separability at the case of C = De~'^^ . Because if a state corresponding to C = _De~^* 
is PPT separable, then another state with stronger phase damping (with increasing 7 while preserving all other 
parameters unchanged) is definitely PPT separable, for this stronger phase damping state we have C > De~^* 
and it is PPT separable. Our proof of the PPT separability of the state at C = De~^* is not a most general 
proof. Wc can only prove the non-negativity of M*^™) by algebraic programming up to m = 17 at the case of 
C = De-^\ Denote m/™^ = {A^ - B'^)-^C"'nI^\ and let A^t^J) (with its elements , < < j) be 

the j-th main submatrix of N^"^\ then to prove the non-negativity of M^™) is to prove that the determinants 
of aU Ar(™j) are not negative. The algebraic programming gives det A^(™^^^ = d-P{d - l)J(J+i)/2p(™J) (d) > 0, 
where d = D/C = 6"^* > 1, and P('"'J)(d) is a polynomial of d with all its coefficients being positive integer, p is 
some integer rely on j. The algebraic programming runs for all m < 17 and proves the non-negativity of M^™). 
We suggest that M^™'' is also non-negative for m > 17, but this is not verified because of the computing time 
of the algebraic programming. 

Another direct way of proving comes from perturbation theory. Firstly A^f^™^ can always be symmetrized. 
The zero order matrix is M^^'^ = M''"^(7t = 0) which is just the case of gaussian state, and all its eigenvalues 
and eigenvectors are well known. So that the first order and second order perturbation of the eigenvalues of 
can be obtained. A more concise way to obtain the perturbation result is as follows: In the eigenequation 

of characteristic function, if we use {a\ $) = exp(— i \a\'^)J2n=o'^'^"^\'^i + 02)™ "(o^i ~ 0^2)" ^ test wave 
function, we then get a matrix (see Appendix) 



/-»pi-$iE(7::)(:) <») 



{C + Dcosx)"'-"-^+''{C - Z3cosa;)'=(-l)'-'^(iDsina;^'+"-2fc 



which is a linear transformation of Af'™^ and has the same eigenvalues. The zero order of Af^™^ is the matrix 
M[^™^ = Af'(™'(7i = 0) which is a diagonal matrix. Hence eigenvalues up to the first order perturbation of 

]\^'{m) gjj^piy Mn^\ When is quite small, Mnn^ can be approximated as Mn^^ ~ {A"^ — B'^)~^{C + 
D)"'-"{C - D)"f{m, n) with 

D D 2_D2 

f{m, n) = 1 - lt[^-^{m - n) + ^-^n + ^ c){D - C) ^"^ " "^'l- ^^^^ 

For odd n, the n — th eigenvalue of the zero order approximation is negative. If f{m, n) is also negative, then 
the n ~ th eigenvalue of the first order approximation becomes positive. For given and D/C we can always 
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A-l-Bexp(-yt) 
- • C-Dexp{-yt) 
- - A-1-B 
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t 

Figure 1: The upper solid line is for logarithmic negativity (LN), the down solid line is for coherent information 
(CI). The parameters are r = 7 = n = r = 0.5. 



find sufficient large m and proper n so that f{m,n) is negative, hence the eigenvalues with large m become 
positive faster than that with small m under phase damping. We need to know at what condition all original 
negative eigenvalues become positive or zero. It is easy to obtained that when m = rt = 1, f(m,n) reaches its 
maximum. Hence when /(1, 1) = 1— ^tjj^ < 0, that is 



C>Dil-jt), 



(11) 



all other odd n negative eigenvalue become positive. Inea. Hll| ) is the first order approximation of Ineq.®. 
Hence at the sense of first order approximation Ineq.® is sufficient for a state to be PPT separable. 

As a by product, we can use matrices M^™"^ to calculate the logarithmic negativity which is an entanglement 
measure itself and provides an upper bound to the distillable entanglement It can be seen in the figure 
when t > t2,the logarithmic negativity is zero, while it is positive when t < t2. For sufficiently small 7t, the 
negativity of the state can be estimated. It is the absolute of the summation of all eigenvalues with odd n. The 
result will be J\f(ti\ w - -yf D(i-C+D) l-A+B _ . B 

lesUlL wm Ui^JW^pj ~ i+C-D l'' {1-C-D)(l+C-D)^ ~ 2(A-B)-l '''[2(A-B)-1]2 • 

One of the most important quantities of a state p is its entropy S [p) — —Trp log p. The entropy of our damped 
state can be obtained by solving the characteristic function eigenequation which is / ^^^x(//, H* ,t) {(3\ D{—fi) 
\a) {a\ $) = A {(3\ <&) . After integrals on p,, a and x, then compare the coefficients of each /3* item of the two 
side, one can get a series of matrices L(™) whose eigenvalues are that of the damped state p and 



-(m) 



B2 ^ 

k 



m 

n - 



^2k jjy+n-2k ^-lt(n-l)^ 



(12) 



log(L(™)). The reduced state 
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The entropy of the state p will be S{p) = -Trit") log(L(°)) - 2YZ=iTrL^ 

of p \s pi — Tr2p with its characteristic function xi(/ii,i) — x(/ii,0, i) — exp[— A' |//|^], hence its entropy is 
S{p\) = A\ogA—{A—l)\og{A—l). It has been proven that the coherent information P — max(0, S'(pi) — S'(p)) 
provides lower bound of distillable entanglement of the state jEl- The coherent information is calculated and 
plotted in the figure. At time the coherent information turns to zero. In the figure we have < ti, but for 
other parameters we may have tg > ti. 

We have investigated the symmetric damping setting of the two mode squeezed thermal state, that is, the 
two modes undergo the same damping and noise. The generalization to asymmetric damping setting and x — p 
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symmetric Gaussian state is straightforward with the method developed here. Denote = 1 — 'J~^Z\b\^j 



D = — —B-^ the M(™) matrix will be 



- . l+n-2k 



The whole issue of the positivity of the asymmetric setting is equivalent that of symmetric setting and omitted 
here. The necessary and sufficient criterion of the PPT separability of the state will be 



VC^>\D\e-^\ (14) 

and the PPT sufficient criterion is again obtained at the sense of algebraic programming and perturbation 
theory. Ineq.® will be generalized to 



- 1){A2 - 1) > \B\ e-^\ (15) 



4 Separability 

We will prove that 

- 1)[A2 - 1) > \B\ ^ > \D\ (16) 

is the sufficient condition of the separability of the damped state p. Let us first consider a Gaussain density 
operator pQ with its characteristic function xg = exp[— ^'j^ Imi|^ ~ ^2 ImsI^ + -B/ii/i2e" + B*/iJ/i2e~")]. The 
Fourier transformation of XGexp(i \p^) is a probability distribution function (pdf) if -^Z [Ai — 1)(A2 — 1) > 
jSjlS', where x is absorbed into /i.This pdf enables the P-representation of pa- Hence pc is separable when 
VC1C2 > \D\. The P-representation of p is a positive integral of the P-representation of pa- Thus p is 
separable. From physical consideration, we may think p is the phase damping of pc, thus when pQ is separable, 
p should be separable. 

The problem left is that when \/CiC2 < \D\ and \/CiC2 > \D\ e^'^*, the state is separable or not. We have 
strong evidence to elucidate that the state is not separable, although we do not give a full proof. The evidence 
is like this: the state can not be expressed in P-representation for \JC\C2 < \D\ , p — J P{ai, 02) |aiQ;2) (aia2| 
cPai(Pa2/TT^ is only possible for VC1C2 > l-D] , where P(ai, 02) is a pdf and \a1a2) denotes two-mode coherent 
state. 

The Fourier transformation will be P{a) — /x(m) 6xp[^ — pa* + p* a]d'^ p / tt'^ . After the integral of p, 
we have 

P{a) I cxp[-E2 |ai|^ - Ei \a2\^ + Fe"aia2 + F*e~"ala*2], (17) 

where c = ^A^_^)f^l-i)-\B\^ = (M-i)ti;-\)-\B(^ ^ ^ = (A,-i)(Af-i)-|Br denote a, = r,e'''%F = |F| e^-^ 
and ip = 9i + 62 + (j), then 

P{a) = cexp[~E2rl - E^rl] £ (^^j exp[-7t |2/ - n|' + i{2l ~ n)^]. (18) 

n=0 ■ i=0 ^ ' 



Clearly P{a) is real, and P{a) is positive. The positivity of P{a) is warranted by the fact that when yJCiC2 > 
\D\ , the state is separable, P{a) is a pdf. If we fix F while decreasing Ei to reach a state with \/CiC2 < \D\ , 
the sign of P{a) will not change by decreasing Ei. Thus P{a) is positive even when \/CiC2 < \D\ . 



What left is the singularity of P{a). We will prove that P{a) is singular if and only if \/CiC2 < \D\ . The 
singularity may appear when \ai\ — Vi ^ 00, P{a) — s- 00. The maximum of P(q;) reaches when (p — 0, thus in the 
following, we set ip = 0. Let g{z) = J2n=o T,?=q (?) exp[-7i \2l ~ then we have g{z) « '£n=o n\J\+-tn 

~ 1+Er=i > l + (;iTT)T = 1 + ^7(^^-1^^)' where we have used DeMoirve-Laplace theorem 

« ^^exp[-(^]. Thus 1 + ^i(e^ - 1 - z) < g{z) < e\ lim^o^iii^ = 1. We arrive at 
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-F'('"ij''2) — > cexp[—E2r1 — Eir2 + 2\F\rir2] < cexp[— 2('\/i?2-Ei — |-F'|)f'if2] when rir2 ^ oo. The non-singularity 
condition of P{a) is simply \fE2E{ > \F\ , which is equivalent to Inea. lfT^ . 

We now compare all three conditions of the separability of p. If \/CiC2 > \D\ , the state p is separable, 
needlessly to say we have > \D\ e-^\ and ^{A-y - 1){A2 - 1) > \B\ e""^*; if y/ChCh < \D\ and > 

\D\ , we have yJ{Ai - 1){A2 - 1) > |S| e-T*,hence Ineq.(fTB)l can be dropped as a necessary condition 
because it is weak than Inea. l|14|l . at this case we do not know the state p is separable or not, we suspect that 
the state is bound entangled; if \/CiC2 <\D \ e~'^*, the state p is entangled. The conditions are expressed with 
the curves A — 1 — Be^"''^, C — De^"''^ and A — B — 1 in the figure for the special case oi Ai — A2, B = B* . 
The zero points of the curves are ti,t2,t3, and ti <t2< t^. For channel without phase damping, the state is a 
gaussian state. All three conditions will be the same, the zero points of the curves will coincide and ^1=^2 = ^3- 

5 Conclusions and Discussions 

In conclusion, the phase damping equation of a state is obtained in the form of characteristic function. It turns 
out to be a usual dissipation equation with respect to the phase angle of the complex variable of the characteristic 
function. The time evolution solution is given for any continuous variable state undergo simultaneous amplitude 
and phase damping and thermal noise. Two of the criteria are given for the amplitude and phase damping of 
a two mode x — p symmetric Gaussian state. One is the sufficient condition of the damped state. The other 
is Peres-Horodecki criterion which is not only necessary but also proved to be PPT sufficient. The proof is 
at the sense of algebraic programming and also perturbation theory. The logarithmic negativity and coherent 
information of the damped state are investigated. 

The evolution of the state is like this: the entanglement of the state (if the state is prepared entangled 
initially) decreases with time, at some time it reaches 0, this time is determined by Peres-Horodecki criterion. 
Then the state may be bound entangled at the next time interval, we proved that the state has not a P- 
representation at this time interval. The end of this time interval is the time determined by the sufficient 
condition of the separability. After this time the state is separable. 

For a channel without phase damping, the state remains a Gaussian state. The two criteria will coincide 
For pure phase damping channel, F = 0, n = 0, we can see that the initially two-mode vacuum state will never 
evolve to a separable state. Hiroshima mentioned this result with numerical calculation|S]. 

This work was supported by the National Natural Science Foundation of China (under Grant No. 10347119), 
Zhejiang Province Natural Science Foundation (under Grant No. R104265) and AQSIQ of China (under Grant 
No. 2004QK38) 



6 Appendix 

From the eigenequation of partial transposed density matrix p^^{t) |<I>) = A |$) , the eigenequation for charac- 
teristic function of it can be deduced as / ^^X^^(m, M*, *) {P\ D{-n) \a) {a\ $) = A(/3| $) . Let (a| $) = 



exp(-i \a\^)Y^ 



m (m) *m—n 
n=0 "1 



, then 



mn 



exp(i|/3|2) J ^^x^^(^,/,t) {P\D{-p)\a) exp(-i 





—^exp[-~Ai \pif ~ A2 \p2\^ - S^i/ize" - B>^/Z2e " + pP*] 



*[Pl - /^l) (P2 - /^2) • 



Where the integral formula 
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is used to integrate a. This formula can further be used to integrate fi. After the integral of /ii we have 
Iran - ^/ f ^ A2 \^^2\' - ^ B* f,,e-^^(il - B i^^e^ + (^2^*2] 

Al JxJ Al 

After the integral of [I2 we have 

with = 1 — -. — I and D = —. — . ^ , ^.a . For each m we expand the binomial and complete the integral of 

A\A^--\B\ A\A'2 — \B\ ^ i 

x-i Iran wiU bc a polynomial of and /32- The eigenequation will be 

m m 
n=0 fc=0 

By comparing the power of /3* , and absorbing the phase of D into the coefficient cl™' , we at last get a matrix 
as in Eq.l^nj ) whose eigenvalues are that of p^^ (t). Eq.Q can be deduced in the same way. 
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